arXiv:hep-th/9502096vl 15 Feb 1995 


SUNY-NTG-95/4 


Effective Lagrangians and Chiral Random Matrix Theory 


M.A. Halasz and J.J.M. Verbaarschot 

Department of Physics 
SUNY, Stony Brook, New York 11794 


Abstract 

Recently, sum rules were derived for the inverse eigenvalues of the Dirac operator. 
They were obtained in two different ways: i) starting from the low-energy effective La- 
grangian and ii) starting from a random matrix theory with the symmetries of the Dirac 
operator. This suggests that the effective theory can be obtained directly from the ran¬ 
dom matrix theory. Previously, this was shown for three or more colors with fundamental 
fermions. In this paper we construct the effective theory from a random matrix theory 
for two colors in the fundamental representation and for an arbitrary number of colors 
in the adjoint representation. We construct a fermionic partition function for Majorana 
fermions in Euclidean space time. Their reality condition is formulated in terms of com¬ 
plex conjugation of the second kind. 
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1 Introduction 


Although it has been widely accepted that QCD is the correct theory for strong interac¬ 
tions, the nonlinearities in the interactions have made it very difficult to obtain accurate 
results that can be compared to experiment. In order to obtain rigorous results many 
researchers have studied a parameter range for which the theory simplifies but neverthe¬ 
less contains the essential features of QCD. Well-known examples are e.g., the large W 
expansion ia, the small volume expansion |^, |^, 2d theory 0, etc.. Such methods have 
provided us with numerous insights in the physics of the full theory justifying any new 
proposal in this direction. 

Recently, Leutwyler and Smilga |Q proposed to focus on the quark mass dependence of 
the Euclidean QCD partition function close the the chiral limit (quark masses m <C Aqcd) 
in volumes with length scale L much larger than a typical hadronic length scale (Aq^^) 
but still much smaller than the pion Compton wave length ( l/^mkQCD)- The first 
condition assures that only the low-lying excitations, in particular the Goldstone modes, 
contribute to the partition function, whereas the second condition allows us to ignore the 
kinetic terms in the Lagrangian. The mass dependence of the partition function is then 
given by the Lagrangian 

C = mI/SRe(TrU), (1.1) 

where S is the vacuum expectation value of V’V' which is assumed to be nonzero, and the 
unitary matrix U parameterizes the Goldstone fields. The question that was asked in Q] is 
to what extent the knowledge of the finite volume partition function puts constraints of the 
spectrum of the Euclidean Dirac operator. By expanding both the QGD partition function 
and the partition function corresponding to (1.1) in powers of m, it was found, by equating 
the coefficients, that the inverse powers of the eigenvalues satisfy the Leutwyler-Smilga 
sum rules. Because of the Banks-Gasher formula the smallest nonzero eigenvalue is of 
order 1 /Aq(^j)L^, whereas in the absence of interactions, the smallest eigenvalue is of order 
1/L. The sum rules hold for eigenvalues well below the latter scale, which for sufficiently 
large volumes is well separated from the scale of the smallest eigenvalues. 
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In previous works i, i 0, 0 we have investigated several questions related to the 
nature of these sum rules. In particular, it was argued that since the effective theory 
is solely based on the symmetries of QCD, the sum rules depend on the symmetries of 
the Dirac operator as well. So when one starts out with a theory with the same global 
symmetries as QCD but no other dynamical input one should arrive at exactly the same 
sum rules. We have constructed such theories: chiral random matrix theories, which apart 
from the chiral symmetry and possible anti-unitary symmetries also contain a remnant of 
the topological structure of QCD. In the framework of random matrix theory one has three 
different universality classes, those with real, complex or quaternion real matrix elements 


T^. For SU{2) with fundamental fermions the matrix elements are real, they are complex 


for more than two colors with fundamental fermions, whereas they are quaternion real 
for adjoint fermions. Indeed, it was shown that the sum rules that are obtained from 
the chiral random matrix theories [|l^ coincide with those obtained from the effective 
theory i0. It should be noted that the sum rules for SU (2) were found only after the 


introduction of chiral random matrix theories [jl^- The three classes of random matrix 
theories correspond to the three different schemes of chiral symmetry breaking, which 


were discussed before in the literature [13 


This suggests that it is possible to derive the hnite volume partition function directly 
from the random matrix theory. In fact, this task has been performed for QCD with three 
or more colors |^. However, for QCD with two colors or QCD with adjoint fermions, the 
situation is more complicated, and up to now such equivalence has not been proved. The 
main objective of this paper is to show that also in these two cases there is a one to one 
correspondence between chiral random matrix theory and the low energy hnite volume 
partition function. 


The main complication is the presence of anti-unitary symmetries in the QCD La- 
grangian. For = 2 this symmetry leads to a real Dirac operator and apart from 
a somewhat more complicated algebra it is straightforward to obtain the effective La- 
grangian. To deal with adjoint fermions we have to face the well-known assertion that 
Majorana fermions do not exist in Euclidean space time [H|. However, this statement 
refers to the transformation properties of a Dirac spinor under Lorentz transformations. It 
does not exclude the possibility to write down a partition function for Majorana fermions 
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in terms of Grassmann integrals (see, for example, |T^). Starting from the observation 
that the adjoint Euclidean Dirac operator is anti-symmetric up to a charge conjugation 
matrix, we indeed succeeded to do this. In fact, using conjugation of the second kind 


Tq, pG|, pT^, the Majorana condition in Euclidean space time is completely analogous to 


the one in Minkowski space time. 


In a previous work |j^, it was observed that the simplest sum rule for each of the 
three cases could be summarized by one formula involving the dimension of the Goldstone 
manifold of the theory. In order to show that this was no coincidence, we present a 
derivation that leads to this result naturally. 

The structure of this paper is as follows. The symmetries of the Dirac operator for 
fundamental and adjoint fermions are discussed in sections 2 and 3, respectively. In 
section 4 we discuss the random matrix theory with these symmetries as input. The 
effective theories for the three different cases are derived in sections 5a, 5b and 5c. A 
general derivation of the simplest sum rule is given in section 6 and concluding remarks 
are made in section 7. 


2 Symmetries of the Dirac operator for fundamental 
fermions 


In this section we study the symmetries of the Euclidean Dirac operatoiQ 


i'lf.Df, = ( 2 . 1 ) 

for a hxed background held in the fundamental representation of SU{Nc). For > 3, 
the Nc X Nc matrix is complex valued, and the only symmetry of the Dirac operator 
is the chiral symmetry 

_ = ( 2 . 2 ) 

^Our conventions are that the Euclidean gamma matrices are Hermitean and satisfy the anti¬ 
commutation relations We use a chiral representation in which 75 is diagonal. 
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Because of this the eigenvalues occur in pairs ±A or are zero. The eigenfunctions are 
related by (j)\ = 750_a- For A = 0 we have the possibility that 750 a=o = ±0a=O) resulting 
in a fermionic zero mode that is either right handed or left handed. It is well known that 
this happens in the held of an instanton [|T^ . 


Because of the symmetry (2.2), the matrix representation of the Dirac operator sim- 
plihes in a chiral basis 0Lfc,} with = (t>Rk and 750Lfc = -4>Lk- The action 

becomes 

\T ^ y r \ :k / 

Dr 


Nf 

EE 

/,g=l kl 


f 

Xr 

f 

Xl 


J^RL 
Dlr imfg 


kl 


Xr 

Xl 


(2.3) 


where we have also included a mass matrix m*(l + 75)/2 + m(l — 
The matrix elements of given by 


/ d^x(j)*^khD(j)Li (2.4) 

are havor diagonal. For a Hermitean Dirac operator we have Dm = In fhe presence 

of fermionic zero modes the number of left handed and right handed states is not necessary 
equal. In general, the matrix Dji^ is a rectangular one. Indeed, the total number of zero 
eigenvalues of the matrix in (2.3) in the chiral limit (m —> 0) is given by the absolute 
value of the difference in dimensionality of the left handed and the right handed Hilbert 
spaces. 

For Nc = 2 the Euclidean Dirac operator is given by 

ilr^Df, = i'jf.dg + (2.5) 

where are the Pauli spin matrices {tiT 2 = etc.). In addition to the chiral symmetry 
(2.2), this Dirac operator possesses the anti-unitary symmetry: 


[zyD, Ct 2 K] = 0 for = 2. ( 2 . 6 ) 

Here, C is the charge conjugation matrix (C = 7274 and = —1) and K is the complex 
conjugation operator. The anti-unitary operator CT 2 K satishes 

{CT2Kf = 1 (2.7) 


From the analogy with the time reversal symmetry in quantum mechanics it is clear 
that this condition allows us to choose a basis in which i'jD is real. To show this, we 
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construct a basis {'ipk) that diagonalizes CT 2 K with eigenvalues one, 


Ct2'iPI = iJk- 


( 2 . 8 ) 


We start with an arbitrary basis vector 0i, then -01 = 0i + Cr20i satishes ( |2.8|) . Next 


we choose 02 perpendicular to 0i. The second basis vector satisfying (p.8|) is given by 
02 = 02 + C'r202) cind it can be shown easily that it is orthogonal to 0i. The next basis 
vector is obtained from 03 perpendicular to both 0 i and 02 , etc.. 

By using (2.6) in the form 


T2Ci'yDCT2 = -{i'jD)*, 


(2.9) 


for which no analogy exists for N^. > 3, it follows immediately that the matrix elements 


(2.4) of the Dirac operator are real in this basis. Note that, because [ 75 ,Cr 2 i^] = 0, the 
above construction holds true for a chiral basis. For Nc = 2, the fermionic part of the 
action is therefore given by (2.3), but with real valued matrices and D^ji. This 
makes it possible to rewrite (2.3) for m = 0 as 



( 2 . 10 ) 


showing that the chiral symmetry group is enlarged to U{2Nf). A mass matrix propor¬ 


tional to the identity breaks this symmetry to Sp{2Nf). Below we will show that the 


same breaking pattern occurs by the formation of a quark condensate. It is important 


that the construction of the above basis relies on the symmetries of the Dirac operator 


only, and that the Dirac matrix becomes real for an arbitrary SU ( 2 ) color gauge held. 

3 Symmetries of the Dirac operator for Majorana 
fermions 

The situation with Majorana fermions in the adjoint representation of the gauge group is 
much more complicated. The reality condition that dehnes Majorana fermions and their 
partition function in Minkowski space cannot be simply generalized to Euclidean space 
time, and in the literature one can End statements that Majorana fermions do not exist 
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in Euclidean space [Q . Below we will make this more explicit and provide a construction 
of a fermionic partition function that works as well in Euclidean space as in Minkowski 


space. 

The Euclidean Dirac operator in the adjoint representation given by 

ilf^D^ = i'yf,{df,6ab +fabcA^), (3.1) 

where is an SU (iVc) background held distributed according to the gluonic action. The 
essential difference from the Dirac operator in the fundamental representation is that the 
long derivative is anti-symmetric under transposition, which is also true in Minkowski 
space. Below we will exploit this by using that the full Dirac operator is anti-symmetric 
up to a constant matrix with determinant one. 

Apart from the chiral symmetry, 

{*7m^m,75} = 0, (3.2) 

the Dirac operator ( p.l|) also has an anti-unitary symmetry 

[i'yD,CK]=0. (3.3) 

First notice that, because 

{CKf = -1, (3.4) 

it is not possible to repeat the construction for A), = 2 in the fundamental representation 
resulting in a real matrix. To see this, remember that the construction was based on the 
diagonalization of the anti-unitary symmetry operator. Now, let us assume that we can 
hnd an eigenvector CKcj) = \(f). Then, 

{CKfcj) = CKX(j) = (3.5) 

which in view of ( p.4|) leads to an obvious contradiction. As a corollary it follows that 0 
and CK(j) are linearly independent. 

To analyze the implications of the symmetry ( p.3|) in Euclidean space, we hrst address 
the issue of obtaining a fermionic action for Majorana fermions in Minkowski space. In this 
case the Dirac operator in the adjoint representation satishes the commutation relation 

[i7^D^,72A'] = 0, (3.6) 
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and the anti-unitary charge conjugation operator 


(72A-)" = 1. 


(3.7) 


Therefore, the operator 72 ^^ can be diagonalized, and a superselection rule can be im¬ 
posed that restricts the partition function to states with eigenvalue 1. Such states, called 
Majorana fermions, can be parameterized by 


ij^ = 


Xr 

-(^2X*R 


(3.8) 


Since 72 ^^ also commutes with Lorentz transformations^], transforms as a Dirac spinor 
under the Lorentz group. This result is usually stated as follows [j^, xl and —(T 2 X*r 
transform in the same way under Lorentz transformations. The Majorana Lagrangian is 
thus given by 


, (3.9) 

where = 74 '^^*. Consistent with 
basis (|3.8|) . 

Let us pursue the construction of Majorana fermions and their action along a different 
route. The crucial observation is that the matrix 7274 * 7 ^ is antisymmetric under trans¬ 
position. This allows us to write down a fermionic action with half as many degrees of 
freedom 

det^^^fyD = det^'^^ 7274 * 7 /^ = J V-ip exp['ip'y 2 'y 4 'i'yDip]. (3.10) 

Lorentz invariance follows immediately form the fact that that 7274 '*/' and 74 '?/’* transform 
in the same way under Lorentz transformations. This can be made more explicit by 
introducing a conjugation operation such that 


(3.7), the Minkowski Dirac operator is real in the 


74 V'* = 7274^, 


(3.11) 


which, of course, is precisely the Majorana condition. Before attacking the problem 
of 4 dimensional Euclidean Majorana fermions let us hrst consider the problem in 1-1-1 
dimensional Minkowski space and two dimensional Euclidean space. 

^In a chiral basis they are given by Q ipL —> X^ipL and ipR —> KriPr with = exp(|(T(a; — iv)) 
and Kfi = (72 1 x 2 . 









In 1+1 dimensional Minkowski space time, the Dirac operator in the adjoint repre¬ 
sentation, i'y^D^ with 7 —matrices dehned by { 7 ^, (We nse the convention 

7 o = <^ 2 , 7 i = icTi) has the anti-nnitary symmetry 

[i'^D,K]={). (3.12) 

Since K also commntes with Lorentz transformations^ the Lagrangian for Majorana 
fermions, dehned by the condition = - 0 ^, is given by P 2 [] 

= (3.13) 

Here, ^jJ^ = = 7o'0^- 

As in 3+1 dimensions, we could have followed an alternative route leading to the same 
Majorana action. The starting point is the observation that ia 2 i'jD is antisymmetric 
with respect to transposition which allows us to write the square root of the fermion 
determinant as a Grassmann integral with only half as many degrees of freedom: 

det^^^iyD = det^^^icr2i7D = J Vijj exp[ipia2i'yDip]. (3-14) 

This expression is Lorentz invariant. Indeed, '0 satishes the Majorana condition. 

Let us now proceed to Euclidean Majorana fermions in two dimensions. The gamma 
matrices are dehned by { 7 ^, 71 /} = and we use the representation 70 = (T 2 , 71 = cxi. 
The Dirac operator satishes the commutation relation 

[i'y^Df,,ia2K] =0. (3.15) 

Because {ia 2 KY = —1 we cannot impose a Majorana condition of the form ia 2 K'ijj = 'ip. 
However, the matrix ia 2 i'jD is anti-symmetric with allows us to halve the number of 
fermionic degrees of freedom in the partition function 

det^^^iyD = det^'^^icr 2 f 7 D = J Dip exp[ijjia2i'yDijj]. (3.16) 

This action is invariant under Euclidean Lorentz transformation, i.e. under ip —^ exp{ia 3 (p)ip. 
This can be made more explicit by introducing a conjugation operator such that 

Ip* = —ia2ip. (3-17) 

^For our representation of the gamma matrices, Lorentz transformations of spinors are given by ^ ^ 
exp{-Lua3)ip. 
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In components, the equations read 


■lAi" = -ip2, i'i = ■i/'i, 


(3.18) 


which forces us to impose the consistency condition 'ipl* = —'ijjk- This conjugation, called 
conjugation of the second kind, is well-known in the mathematical literature on Grass- 


mann variables [jTg , and has been used extensively in the supersymmetric formulation of 
random matrix theories |]^ |^]. The Majorana constraint (|3.17|) 

V'l \ f '4’i\ 


i(J2K 


(3.19) 


\ V'2 J V '^’2 I ' 

makes the symmetry ( 11.15 ) manifest at the level of the partition function. 

Let us now return to 4 Euclidean dimensions. The strategy should be clear: we 
construct an anti-symmetric operator that differs from the Euclidean Dirac operator only 
by a factor with unit determinant. As can be seen from the following theorem, it is not 
an accident that this works. 

Theorem. Consider a Hermitean operator H such that [if, AK] = 0 and A = AA^ = 1. 
If {AKY = —K^ then (HA)'^ = —HA. If, moreover A or iA are orthogonal then also 
{AHf = -AH. 

The proof of this theorem is immediate. It can be applied to the Dirac operator 
in Euclidean space time. However, since going from Euclidean to Minkowski space just 
amounts to multiplying the spacial gamma matrices by a factor i, in both cases the Dirac 
operator behaves the same under transposition. 

Using this theorem, one concludes from (3.3) that Ci'yD is an antisymmetric matrix, 
which allows us to write down a fermionic partition function with only half as many 
degrees of freedom 


det^'^^iyii = det^'^^C'iyii = J 'D%jjexp[ipCi'jD%jj]. 


(3.20) 


This construction also works in the presence of a mass term. It is straightforward to 
verify that (3.20) is invariant under Euclidean Lorentz transformations. As before, this 
can be made more explicit by choosing the components of ip such that they satisfy the 
conjugation equation 


Cip* = ±ip. 


(3.21) 
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A solution with eigenvalue +1 can be parametrized by 


^ = 


( \ 

'lp2 

V -^2 / 


(3.22) 


which can be viewed as a Majorana constraint on the fermionic integration variables. As 
was the case for Euclidean fermions in 2 dimensions, we have to impose conjugation of 
the second kind on the Grassmann variables (otherwise the equation CK'ijj = ±-0 does 
not have solutions, see eq. (3.5)). 

To display the matrix structure of the Dirac operator we express the held in the La- 
grangian (|3.20|) in terms of a complete set with Grassmannian coefficients. An expansion 
consistent with the Majorana constraint (|3.21|) is given by (because {CKY = —1, the 
c-number functions 0 and C(f)* are linearly independent) 


(3.23) 


(3.24) 


'4’R-YI (t>RkXRk + 
k 

for right handed fermions and, 

'ipL=Yl (pLkXLk + C(t)lkX*Lk 

k 

for left handed fermions. In this basis the fermionic piece of the action is given by 

f Xr ^ 

xh 

/ 

Xl 

f * I 

\Xl J 

where each 2x2 block of Drl is given by 

(r, X I(p*Rk'^xD(f)Li \ . . 

{VRL)ki jct>nkC^tlDct>Li j4>RkC^tlDCct>li ) ’ 



( \ 


q * 

Xr 

^LR ^fg 

Xl 


\ xV / 


(3.25) 


and Dlr = Drl- (|3.25|) we have included the havor indices. Using that 

C^ijDC = {ixDy, 

C^ijD = -{ijoyc, 


(3.27) 

(3.28) 


we hnd that each 2x2 block in is quaternion real, i.e. of the form aQ + iakCk with 
real. The mass matrix is proportional to the unit quaternion but is generally not diagonal 
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in the flavor indices. Since the product of Grassmann variables that multiplies the mass 
matrix is symmetric in the flavor indices, the mass matrix can be taken symmetric as 

well. 


Because the Grassmannian vectors on the left and on the right in (3.25) do not trans¬ 
form independently under transformations in flavor space, the chiral symmetry for m = 0 
is reduced to U{Nf). For a nonzero mass matrix that is a multiple of the identity, only an 
invariance under the 0{Nf) subgroup remains. Below we will see that the same breaking 
is achieved by the formation of a nonzero chiral condensate. 


The eigenvalues of the Dirac operator (|3.25|) are doubly degenerate and, for zero mass, 
occur in pairs ±A or are zero. The number of zero eigenvalues is equal to twice the 
absolute value of the difference of the dimensionality of the right handed and the left 
handed space. 


4 Chiral random matrix theories 


The matrix elements of the Dirac operator fluctuate over the ensemble of gauge fields. In 
this section we introduce a model partition function with the symmetries of the Dirac 
operator but with matrix elements fluctuating according to independent gaussian distri¬ 
butions. As discussed in the introduction, we expect that certain low energy quantities, 
such as the fluctuations of the eigenvalues on a microscopic scale, do not depend on the 
dynamics of QGD interaction and can be calculated in this simplified model. 

The QGD partition function can be written as 

Zqcd = ZQCDitr), (4.1) 

V 

where the partition function in the sector with v fermionic zero modes is defined by 

Nf 

Zqcd{^) = {Y[m} n (^n + \mf\‘^))A- (4.2) 

/=1 A„>0 

The eigenvalues of the mass matrix are denoted by m/, and the average goes over all gauge 
fields with topological charge u weighted according to the gluonic action. For fundamental 
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fermions v = but for adjoint fermions, with each of the doubly degenerate eigenvalues 
included only once in the fermion determinant, we have V = NcV (see Q). 

In sectors with a nonzero total topological charge, the number of right handed modes 
in (|2.3|) and (|3.25|) differs from the number of left handed modes. The random matrix 
partition function with n and n + P of such modes, respectively, and Nf flavors is dehned 
by 


Nt 


Zp{u,Nf)= VTll det 

J /=! 


m*f iT 


iTt 


ruf 


exp[ 




TrTT'^], 


(4.3) 


where T is an u x (n + P) matrix. The integration over T is according to the Haar 
measure. The matrix elements of T are real for /3 = 1 corresponding to QCD with Nc = 2 
in the fundamental representation. They are complex for iV). > 3 in the fundamental 
representation (/3 = 2). For gauge helds in the adjoint representation the matrix elements 
of the Dirac operator are quaternion real (see (3.27)), and the matrix elements Tjj in 


( [4.3|) are chosen quaternion real as well (/3 = 4). If we write the quaternions in terms of 
2x2 matrices, the matrix T is a 2u x 2(u + P) matrix. Of course, the masses ruf are 
multiplied by the quaternion unit matrix. The determinant in (|4.3|) for /? = 4 is the so 
called Qdet and the trace is the QTr. It can be shown that for a quaternion real matrix 
A that Qdet^A = det A [^. In a 2 X 2 matrix representation of the quaternions the 
QTr is just one half the ordinary trace. Note that C times the unit in the flavor indices 
multiplied by the matrix in (3.26) is antisymmetric which makes the square root of its 
determinant is well dehned. 

In analogy with the classical random matrix ensembles, these ensembles will be called, 
the chiral orthogonal ensemble (chGOE), the chiral unitary ensemble (chGUE) and the 
chiral symplectic ensemble (chGSE), for /? = 1, 2, and 4, respectively. 

We will identify the total number of modes, 2n, (we always have P -C u) with the 
volume of space time. This corresponds to choosing units in which the density of the 
low-lying modes is equal to one. Below we will see that the parameter S can be identihed 
as the chiral condensate. 

The matrix ensembles described by the partition function ( |4.3| ) are equivalent to what 
is known in the random matrix literature as the Laguerre ensembles. They hrst were in¬ 
troduced by Fox and Kahn [^]. The simplest case {3 = 2 was analyzed in |^. An analysis 
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of the /3 = 1 case in the context of the microscopic spectral density of the Dirac operator 
was given in An analysis which also includes many other correlation functions for 
13 = 1 and (3 = A was performed in a series of papers by Nagao and coworkers and For¬ 


rester 


1. Other results for the chiral random matrix ensembles have been obtained 


in terms of a supersymmetric formulation |27| The recent work on this subject is based 


on results by Mehta and Mahoux who introduced the skew orthogonal polynomials 
originally invented by Dyson p^ . 

The fermion determinant in (4.1) can be written as a Grassmann integral 

Nf 


Z0{V,Nj) = J VTVrVijexp[-t 


f,9=l 


T 


V'f - ^^^^TrTT^]. (4.4) 


ik 


For adjoint fermions {f3 = 4) the matrix elements of T are quaternion real, and each 
component 'ijj{ is a vector of length four given by (see (3.26)) 



(xL \ 

X*r{ 

f 

xii 


(4.5) 


* / 

V XL i / 

They satisfy a reality condition similar to the Majorana constraint. In order to assure 
ourselves of a positive dehnite fermion determinant for i/ = 0 we have included a factor 
i in front of the fermionic action. For (3 = 1 and (3 = 2 the components of and 'ip* are 
independent integration variables. 


5 Effective theory 


In this section we derive the effective theory corresponding to the random matrix partition 
function ( |4.4|) . We proceed by averaging over the matrix T resulting in a four-fermion 
interaction which can be made gaussian at the expense of a new bosonic integration 
variable. After performing the Grassmann integrals the resulting theory is amenable to a 
saddle point approximation in which the integrals over the soft modes are kept and the 
integrals over the hard modes are done to gaussian order. 

For completeness we start the discussion with the simplest case (3 = 2 (section 5a) 
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which was already analyzed in ref. 
and /9 = 4, respectively. 


lOj. In sections 5b and 5c we discnss the cases /3 


1 


5.1 Effective theory for /3 = 2 


After averaging over the matrix elements of the Dirac operator the partition fnnction 
becomes 




(5.1) 


Here and below we have nsed the ~ sign in order to indicate that constant factors have 
been absorbed in the normalization of the partition fnnction. The term of fonrth order 
in the Grassmann variables can be rewritten as the difference of two sqnares, 


i'h'lik'H'ii'Ri = + i’Rii’Rif - - V’KV'D' 


(5.2) 


Using the Hnbbard-Stratonovitch transformation, each of the sqnares can be linearized 
by introducing an additional Gaussian integral 


exp{—AQ^) ~ J (icrexp(—— — iQa). 


(5.3) 


This results in the partition function 


Z2{v-, Nf) 


VaiDa 2 Vil)Vil)* exp[ 




Tr(aicrf + cr2aj) 




r/S 


^fg 


f*. 


^fg 


Jg 


where cxi and a 2 are arbitrary real NfXNf matrices. The integration over the Grassmann 
variables yields 

Z 2 (h, Nf) ~ J VAdet'^~^‘'{A^ + m*)det”(A + m) exp[-^-^TrAA'^], (5.5) 

where A is an arbitrary complex matrix. It can be diagonalized according to 


A=UAV-\ 


(5.6) 


In order to have the same number of degrees of freedom on both sides of the equation we 
chose [/ G U{Nf) and V G U{Nf)/{U{l))^f. All matrix elements of the diagonal matrix 
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A are real non-negative. For n —>■ cx) and ||ni||S -C 1, the integral over the eigenvalues 
can be performed by a saddle point approximation at m = 0. The solutions of the saddle 
point equation for A are given by 

Ai = ±i, (5.7) 

but only the solution with all signs positive is inside the integration manifold. At this 
point, the integral only depends on the combination UV~^ which allows us to absorb V in 
the integration over U and perform the V integration. For small masses, m/S -C 1 (the 
eigenvalues of the mass matrix are denoted hj rrif), the integral over U is soft and cannot 
be done by a saddle point method. However, in this limit we can expand the determinant 
to hrst order in m which leads to the partition function 

Z2(P,Nf) ^ [ deexpUpe) [ T>Hexp[nSTr(m*He*®/^^(5.8) 

J Juesu{Nf) 

where we have split the integration over 17 in a 17(1) integral over 6 and an integral over 
SU{Nf) [U ^ U exp[z6'/A^/]). For a diagonal mass matrix the condensate is given by 

iWf) = ^dmf'^ogZ, (5.9) 

where the differentiation is with respect to one of the quark masses. The right hand side 
should be evaluated for uifuZ S> 1, which allows us to use a saddle point integration for 
U with the result that 

(g/g/) = Ecos(^). (5.10) 

This completes the calculation of the partition function that, with the identihcation of 2n 

as the volume of space-time, was the starting point of ref. p. 


5.2 Effective theory for (3=1 


For /? = 1, the overlap matrix T is real. After averaging over T the partition function is 
given by 


Zi{p,Nf) ~ J 


exp 


1 

2nS2/? 


ijRi 



m. 






'4’Lk \ 

i^Lk If, 
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(5.11) 





The introduction of anti-symmetric unit matrix 

I = 


0 -1 

1 0 


(5.12) 


allows us to rearrange the fermions in multiplets of length 2Nf, in agreement with the 
well-known result that for two colors the the chiral symmetry group is enlarged from 
U{Nf) X U{Nf) to U{2Nf). Since baryons consist of two quarks for Nc = 2, it is not 
surprising that the four-fermion interaction contains both mesonic and di-quark bilinears. 
It can be rewritten as the difference of 


1 

4 






'ipLk 

i^Lk 


'IpLk 

i^Lk 


2 


If! 


f'f 


(5.13) 


+ ^aa’ 

j \ 'T m / g y-rnn/g, \ 

and a similar expression with the plus sign exchanged by a minus sign. Each of the 
fermionic bilinears is anti-symmetric in the flavor indices. Using the Hubbard-Stratonovitch 
transformation (5.3), the squares can be linearized with the help of an anti-symmetric ma¬ 
trix which yields the following partition function 


Zi{P,Nf) ~ J T>V’*^^^exp [ - 2nJ:‘^f3Tr{AA^) 


V’i? 

rn 

n 


(V + Im-) 


fa 


'4^r 

rn 


lAl'^ + -M 


fa 


rL 


The mass matrix is dehned by 


M = 


0 

m 


-m 

0 


(5.14) 


(5.15) 


and y4 is a general antisymmetric complex 2Nf x 2Nf matrix. Carrying out the fermion 
integrals and rescaling Ahy a factor two, the partition function can be compactly written 


as 


Zi{9, Nf) ~ J VAPr+^{A^ + M*)Pr{A + M) exp[-^^Trklklt], (5.16) 

where Pf denotes the Pfaffian of the matrix (For an even dimensional anti-symmetric 
matrix, the Pfaffian is equal to the square root of its determinant with a dehnite choice 
for its phase |^, |^). The complex antisymmetric matrix A can be brought into a 
standard form as 


A = UAU^, 


(5.17) 
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where ?7 is a unitary matrix and A a real antisymmetric matrix with Ak^k+i = —Ak+i,k = 
Afc, k = 1, - ■ ■, 2Nf — 1 and all other matrix elements zero. By redehning U we can always 


choose all > 0. In ( ^.161 ) we use A and U as new integration variables. Note that 
U G U{2Nf)/{Sp{2))^f so that the total number of degrees of freedom on both sides of 
( b.lTp is the same. We are interested in the limit n ^ oo and ||Ad||E <C 1. Then the 
A integration can be performed by a saddle point integration at Af =0, whereas the 
remaining integrals have to be performed exactly for the actual value of the mass. At the 
saddle point inside the integration manifold, i.e. Afc = 1/S, the integrand does not depend 
on U for A1 = 0. For A4 ^ 0 the U dependence is in the form UIU"'", where I is the 
antisymmetric unit matrix. The integration over the soft modes U is thus parametrized 
by the coset U{2Nf)/Sp{2Nj). For ||A1||S 1 the exponentiated determinants can be 

expanded to hrst order in Ai resulting in the low energy hnite volume partition function 

ZAP,Nf) ^ f dOexpUpe) f explnEReie^^/^fTiUIU^M)], (5.18) 

J JueSU{2Nf)/Sp{2Nf) 

which was used as starting point for the calculation of the Leutwyler-Smilga sum rules in 
0 (with the identihcation of 2n as the volume of space-time). The phase of det U has 
been isolated by the substitution U ^ Uexp[i6/2Nf]. Note that in this way the phase of 
the Pfaffian in (5.16) covers the full complex unit circle for 6 G [0,27r]. The integration 
over the stability subgroup only modihes the partition function by a constant. Therefore, 
the integration in (|5.18|) can be extended to SU{Nf), which facilitates further evaluation 


of the partition function |11 


Also in this case the condensate is given by the logarithmic derivative in (^.9|) evaluated 
for uifTiTi 3> 1 (with m/ an eigenvalue of m). In this limit the U integral can be performed 
by a saddle point method. The saddle point is at UIU"'" = 1 resulting in the identihcation 


Q 

{q,q,) =Ecos(—). 


(5.19) 


5.3 Effective theory for /3 = 4 

In this case the overlap matrix elements are quaternion real 
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where ct+ = (—i, a)^ and the are real. In terms of the variables the partition 
fnnction reads 

Z4(z7, Nf) ~ J Va^Vilj*ViP exp [ - ^^^^Tra^aJ + 

+ + (5-21) 

The components of the fermionic variables are 2 component vectors dehned by 


By using the identity ((cr+)l)^ = —cT 2 cr^cr 2 we hnd that the two fermionic terms in (|5.21|) 
are identical. Averaging over the overlap matrix elements leads to 

(5.23) 


Using the Fierz identity 


+ aP +'yS _ 




^ap^'ys') 


(5.24) 


and the representation (^.22 ) of the spinors, the partition function can be simplihed to 

. 4 


Zi{u,Nf) ~ J V^p*Vi) 


exp 


iXRiXRi + XRiXDixilXLk + XlIXli 


npT? 

™-fg'4’R*'^Ri + ^fgi’lki’Lk]- 


(5.25) 


This four fermion interaction can be written as the difference of two squares which can be 
linearized with the help of the Hubbard-Stratonovitch transformation (|5.3|) by introducing 
the bosonic variables ai and ct 2 . Both ai and a 2 are symmetric real valued Nf x Nf 
matrices. In terms of the new variables the partition function reads 

n/3S2 


Z4i^,N,) 


'D'ip*'DipT)aiT>a 2 exp[ 


-Tr(cT^ Taj] 


X exp[2x;|j*xfji(ai - ia2)fg + '^XL*kXLk{(^i + i(^2)fg + 2{m*fgXRiXh + ^fgXilXLk)]- 

(5.26) 

The Grassmann integrals can be performed trivially resulting in the partition function 
Z 4 (n, Nf) ~ I VAdee^‘'{A^ + m*)det"(A + m) exp[-^^TrAA^], (5.27) 
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where ^4 is a complex symmetric matrix. Such a real matrix can be diagonalized by a 
unitary matrix as follows 

A = UAU^ (5.28) 


with A a real nonnegative diagonal matrix. As before we use A and U as new integration 
variables. The mass plays the role of a small symmetry breaking term (||m||S 1) so 

that the integration over the hard modes, the A, can be performed by a saddle point 
approximation at m = 0. The saddle points are given by 

Afc = ±^, (5.29) 


but only the saddle point with all signs positive is inside the integration manifold. The 
integral over the soft modes has to be taken into account exactly, but the exponentiated 
determinant can be expanded to hrst order in m. Writing U as the product of a phase and 
a special unitary transformation {U —>■ exp{i6/2Nf)U) we arrive at the partition function 


ZAv,Nf) ^ ! 1)17 exp[2nSRe(e*®/^^Trf/f/^m)]. (5.30) 

J JueSU{Nf)/0{Nf) 

Because the integrand depends only on U'^U, the integration over U is effectively over 
the coset SU{Nf)/0{Nf), where the volume of the stability subgroup can be absorbed 
in the normalization. As before, the condensate is given by the logarithmic derivative of 
the partition function for m/nS S> 1 (The eigenvalues of m are denoted by rrif). The 
integral can than be performed by a saddle point method resulting in the identihcation 


Q 

{q/q,) = Ecos(—). 


(5.31) 


With the identihcation of 2n as the volume of space-time, the partition function (|5.30|) 
was the starting point for the derivation of the Leutwyler-Smilga sum rules in 


11 


6 Sum rules for the eigenvalues of the Dirac operator 

As already mentioned in the introduction, the sum of inverse powers of the eigenvalues of 
the Dirac operator satishes sum rules. Since they have been derived both from effective 
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field theory and random matrix theory, we want to restrict ourselves to one feature that 
was not clarihed in earlier work. 


To obtain sum rules we expand both the expectation value of the fermion determinant 
and the hnite volume partition function in powers of (We consider the case of a diagonal 
mass matrix with all masses equal to m). For the fermion determinant we obtain 


Z(m) ^ 1 

^ = l + m]V,y 


( 6 . 1 ) 


A„>0 n 

where, for simplicity, we consider the case with all masses equal. Equating the coefficients 
of in (6.1) and of the expansion of the hnite volume partition function (see below). 


leads to sum rules that can be summarized into one formula valid for each of the 
three cases discussed in section 5: 


— E — = 

y2 \2 


( 6 . 2 ) 


4(|i/| + (dim(coset) + 1)/Af^)' 

The volume of space-time is denoted by V. In this section we give a derivation for ly = 0 
showing that this unifying feature was no accident. 

The coefficient of of the ratio Z(m)/Z(0) for z/ = 0 of hnite volume partition 
functions given in (|5.8|) , (|5.18|) and (|5.30|) involves the calculation of integrals of the form 


C(kl) = [ dVTT{VA)TT*{VA), (6.3) 

JVGG jH 

For z/ = 0 the U{1) integral over 6 can be absorbed in the integration over the coset G/H, 
which now becomes U{2Nf)/Sp{2Nf), U{Nf) and U{Nf)/SO{Nf) for (3 = 1, 2 and 4, 
respectively. The matrix A is the mass matrix which is real anti-symmetric, complex and 
real symmetric in this order. Using the invariance of the measure it follows that in all 
three cases C(44) ~ TrAAh For a mass matrix with equal masses in the standard form, 
a 4 ~ 1 for /? = 2 and (3 = 4 and proportional to the antisymmetric unit matrix I for 


(3 = 1. 


To proceed let us introduce generators of the cosets []TT|. They satisfy the orthogo¬ 
nality relations 


Tr tkti 



(6.4) 


and are real anti-symmetric, Hermitean and real symmetric, for (3 = 1, (3 = 2 and (3 = 4, 
respectively. The total number of generators is denoted by M and the dimension of 
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the matrices U and tk is denoted by d. The generators will be chosen snch that ti is 
proportional to the mass matrix. Becanse of the invariance discnssed in the previons 
paragraph and the normalization ( b.4|) , we hnd that 

c{ti) = at2) = --- = atM), ( 6 . 5 ) 


allowing us to rewrite the integral as 

1 ^ r 

«‘i) = Xf S / , d-VTv{Vh)Tv-‘{Vt,). (6.6) 

By expanding U in its generators {U = J^k'^ktk) and using the orthogonality relation 
( p.4|) and the unitarity of U, the integral becomes trivial. The result is given by 

C(«i) = ^vol(G/ff), (6.7) 


The volume of the coset cancels in the ratio Z{m)/Z{0) of the partition functions. Using 
that the dimensionality d = Nf for /3 = 2 and /3 = 4 and d = 2Nf for /3 = 1, we hnd that 


Z{rn) 

'm 


u=0 


1 + 

4 2M ^ 


( 6 . 8 ) 


The extra factor 2Nf appeared because of the difference in normalization between the 
unit matrix and the generators ( |6.4|) . As already discussed in section 4, we have identihed 
the total number of zero modes, 2n, with the volume V of space-time. 

In (|6.2|) the dimension of the coset does not include the U{1) factor, so in the notation 
of (U) we have 


M = dim(coset) -I- 1. 


(6.9) 


By comparing the coefficients of m? in (|6.1|) and ( |6.8|) we indeed reproduce the sum rule 


7 Conclusions 


Starting from a chiral random matrix theory with the symmetries of the Dirac operator 
in an arbitrary background held, we have derived the hnite volume partition functions 
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that were used as starting point in the derivation of the Leutwyler-Smilga sum rules in 
iin. A derivation of the simplest Leutwyler-Smilga sum rules shows in a natural way 
that they depend on the number of Goldstone bosons per flavor. 


Our results explain the miracle that effective Lagrangians and chiral random matrix 
theory produce the same sum rules. The advantage of chiral random matrix theory is that 
it leads naturally to three classes of sum rules. The absence of a third class in motivated 
the studies in with the result that QCD with two colors in the fundamental represen¬ 
tation constitutes a separate universality class. In the present work we have shown that 
then the Dirac operator is real and that the corresponding random matrix theory (chGOE) 
is equivalent to an effective theory based on the coset SU{2Nf)/Sp{2Nf). For three or 
more colors in the fundamental representation the Dirac operator is complex and the cor¬ 
responding random matrix theory with complex matrix elements (chGUE) is equivalent 
to a hnite volume partition function based on the coset SUniNf) x SUl{N f)/SU{Nf). 
The last universality class is for a Dirac operator with adjoint fermions and two or more 
colors. Then the Dirac operator is quaternion real, and in this paper we have shown 
that the corresponding random matrix theory with quaternion real matrix elements is 
equivalent to a hnite volume partition function based on the coset SU{Nf)/0{Nf). 


Apart from two external parameters, the vacuum angle and the mass matrix, the 
random matrix theory and the corresponding hnite volume partition function depend on 
only one dynamical parameter of QGD; the chiral condensate. The hnite volume partition 
function, and thus the chiral random matrix theory, provides us with the mass dependence 
of the QGD partition function. However, this does not imply that all properties of the 
chiral random matrix model are physical. For example, the average level density, which has 
a semicircular shape, is certainly unphysical. This leads to the question which properties 
of the chiral random matrix theory actually determine the hnite volume partition function. 


In the hnite volume partition function the mass occurs in the combination mV, whereas 
in the QGD partition function or in the random matrix theory, the mass occurs as m/A, 
where A is an eigenvalue of the Dirac operator. This suggests that the ehective theory 
is only sensitive to the distribution of eigenvalues on a scale 1/H. The spectral density 
on this scale, also called the microscopic spectral density, is well dehned. From the 
work on quantum chaos and mesoscopic systems (including nuclei) we know that the 
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distribution of eigenvalues on a microscopic scale is universal and is given by random 
matrix theory. The most impressive piece of work in the present context is by Slevin 


and Nagao who studied the logarithm of the transfer matrix of a mesoscopic system 
in a magnetic held for the Hofstadter model. The random matrix theory of this system 
has the symmetries of the Dirac operator for Nf = 0 and three or more colors in the 
fundamental representation. Indeed, by numerical computations, they found that the 
microscopic spectral density is given by the corresponding chiral random matrix theory. 
This led us to the conjecture that the so called microscopic spectral density is universal. 
However, this does not imply that it can be determined from the hnite volume partition 
function. For example, for one havor, in each of the three cases, Goldstone bosons are 
absent and the hnite volume partition function is the same but the random matrix theory 
is different. This shows that the low-lying spectrum of the Dirac operator cannot be 
derived from the complete set of Leutwyler-Smilga sum rules. 

We conclude with the statement that although the hnite volume partition function does 
not determine the microscopic spectral density, there is ample evidence that if combined 
with universality according to the anti-unitary symmetries of the Dirac operator this leads 
to a unique prediction of the spectral density of the Dirac operator on a scale of no more 
than a hnite number of eigenvalues from zero. 
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